A partial wave analysis of the vr tt 7r+ and tt tt^tt^ systems 
and the search for a J^'^ = 1 + meson 

A. R. Dzierba, R. Mitchell, E. Scott, P. Smith, M. Swat, and S. Teige 
Department of Physics, Indiana University, Bloomington, IN 47405 

A. P. Szczepaniak 

Department of Physics and Nuclear Theory Center, Indiana University, Bloomington, IN 47405 

S. P. Denisov, V. Dorofeev, I. Kachaev, V. Lipaev, A. V. Popov, and D. I. Ryabchikov 
Institute for High Energy Physics, Protvino, Russian Federation 142284 



in 
o 
o 

(N 
> 

o 
O 
vn 

(N 

(N 
> 
00 

O 
O 

in 
o 

• 

^ : 
I 

Oh! 



X 



V. A. Bodyagin and A. Demianov 
Nuclear Physics Institute, Moscow State University, Moscow, Russian Federation 119992 

(Dated: February 7, 2008) 

A partial wave analysis (PWA) of the 7r~7r~7r"'" and 'k~'k'^'k^ systems produced in the reaction 
7r~p — » (37r)~p at 18 GeV/c was carried out using an isobar model assumption. This analysis is 
based on 3.0M ■r~-k'^-r^ events and 2.6M 7r~7r~7r"'" events and shows production of the ai(1260), 
02(1320), 7r2(1670) and a4(2040) resonances. Results of detailed studies of the stability of partial 
wave fits are presented. An earlier analysis of 250K n^ir^n^ events from the same experiment 
showed possible evidence for a J^^ = 1 ^ exotic meson with a mass of ^1.6 GeV/c^ decaying into 
pTT. In this analysis of a higher statistics sample of the (37r)~ system in two charged modes we find 
no evidence of an exotic meson. 

PACS numbers: ll.80.Cr, 13.60.Le, 13.60Rj 
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I. INTRODUCTION 

In this paper we present a partial wave analysis (PWA) 
of a high-statistics sample of events corresponding to the 
production of the (37r)~ system produced in tt~p colli- 
sions in two modes: tt^tt'^tt'' and Tr^Tr^Tr"*". This sam- 
ple size exceeds, by at least an order of magnitude, the 
largest published sample size of Stt events to date. Pre- 
vious Stt analyses led to the discovery and/or determi- 
nation of properties of the ai(1260), a2(1320), 7r2(1670), 
7r(1800) and the 04(2040) resonances [1-8]. 

In 1998, the E852 collaboration reported evidence for 
the 7ri(1600), a J^^ = 1~+ exotic hybrid meson with 
a mass of 1.6 GeV/c^ decaying into pir [9, 10]. That 
analysis was based on 250,000 events of the reaction 
7r"p 'K~TT~TT^p collected in 1994. We report on the 
analysis of additional data collected in 1995 including 
3.0M events of the reaction ir^p — > 7r~7r°7r°p and 2.6M 
events of the reaction 7r"p — > 7r~7r~7r"'"p. 

The identification of exotic mesons requires a PWA to 
extract signals. The Stt system provides a particularly 
attractive venue for such searches since earlier work sug- 
gests a rich spectrum of meson resonances. If an exotic 
meson is produced with relatively small amplitude, its 
interference with nearby well-established resonances is a 
sensitive search tool. 

In almost all of the published amplitude analyses of the 
Stt system, the isobar model was employed - a Stt system 
with a particular J^*^ is produced and decays into a di- 
pion resonance with well-defined quantum numbers and 
a bachelor tt followed by the decay of the di-pion reso- 



nance. This assumption is successful in describing many 
features of the Stt system and is motivated by the ob- 
servation that the di-pion effective mass spectrum shows 
prominent resonance production, for example, the p(770) 
and the /2(1275). The di-pion resonances considered in 
this analysis include the /o(980), p(770), /2(1275) and 
the p3(1690). We also include parametrizations of S- 
wave TTTT scattering. 

In this analysis, we observe the ai(1260), 02(1320), 
7r2(1670) and the 04(2040) in appropriate partial wave 
intensities and in their relative phase differences. We also 
searched for, but find no evidence for, the J^'-^ = 1 ^ 
exotic 7ri(1600) in either the tt~tt^7t^ or tt^tt^tt^ mode. 



In order to extract reliable information from a partial 
wave analysis it is important to establish a procedure for 
determining a sufficient set of partial waves. Failure to in- 
clude important partial waves in the series expansion may 
lead to inconsistent results and erroneous conclusions. In 
this paper we describe our procedure for determining a 
sufhcient wave set. That procedure includes a compari- 
son of moments as calculated from PWA solutions with 
those computed directly from data. We emphasize that 
this analysis is similar to that of references [9, 10] in that 
the same isobar model assumptions are made but the 
final set of partial waves used is different. Both analy- 
ses make the same assumptions about coherence between 
different partial waves. It is possible that relaxing these 
coherence assumptions could lead to different results in 
both analyses. 
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A. General features of the (Stt) system 

The Stt system with non-zero charge has isospin / > 0, 
and, since no flavor exotic mesons have been found, we 
assume 1=1. Since a state with an odd number of 
pions has negative G parity, the relationship G = C(— 1)^ 
imphes positive C parity for the (37r)~ system. 

The simultaneous observation of the two (Stt)^ modes, 
TT^n'^Tr'^ and 7r^7r^7r+, in the same experiment provides 
important cross checks. Consider the production of res- 
onance X and its subsequent decay into (Stt)" via an in- 
termediate di-pion resonance. If the intermediate di-pion 
resonance is an isoscalar, then the yield into 7r~7r~7r"'" 
should be twice that into tt^tt^tt'^. Similarly, if the di- 
pion is an isovector, then there should be equal yields into 
7r~7r~7r"'" and 7r~7r''7r''. Since the 7r~7r~7r^ and Tr^Tr'^Tr" 
modes rely on different elements of the detector, any mis- 
understandings in the acceptance would affect the two 
modes differently and ultimately lead to inconsistent in- 
terpretations of the PWA results. This underscores the 
importance of having these two different modes to vali- 
date the analysis. 



B. Exotic mesons within QCD 

Quantum chromodynamics (QCD) predicts a spec- 
trum of mesons beyond the qq bound states of the con- 
ventional quark model. The spin (J), parity (P) and 
charge conjugation (C) quantum numbers of a qq sys- 
tem are: J = L + S where L is the angular momen- 
tum between the quarks and S is the total quark spin, 
P = (-1)^+^ and C = (-1)^+^. The J^^ combinations: 
0^ , 1 ^, 2^ , . . . are not allowed and are called exotic . 
Lattice QCD [11] and QCD-inspired models predict that 
the excitation of the gluonic field within a meson leads to 
hybrid mesons, where the gluonic degrees of freedom al- 
low hybrid mesons to have exotic quantum numbers. The 
observation of hybrid mesons and measurement of their 
properties provides experimental input necessary for an 
understanding of quark and gluon confinement in QCD. 



C. Organization of this paper 

This paper is organized as follows. Section II is a dis- 
cussion of experimental details, including a description 
of the apparatus (HA), event reconstruction and selec- 
tions (II B) and distributions in effective mass before and 
after acceptance corrections (11 C). An overview of the 
PWA methodology, including selection of the wave set, 
is given in Section III. Section IV is a presentation of 
the PWA results for individual partial waves using the 
wave set used in this analysis and the wave set used in 
the analysis of references [9, 10] for both the 7r~7i-°7r" 
and 7r^7r^7r+ modes. Section V focuses on the exotic 
+ waves and how the failure to include im- 
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FIG. 1: E852 apparatus and the MPS at BNL. 

portant J^'^ =2 waves could lead to misleading con- 
clusions regarding the existence of an exotic resonance. 
Section VI discusses studies of systematic uncertainties 
and the conclusions are presented in Section VII. 



II. EXPERIMENTAL DETAILS 

A. The experimental apparatus 

The E852 experiment used the upgraded multiparticle 
spectrometer (MPS) at the Alternating Gradient Syn- 
chrotron (AGS) at Brookhaven National Lab (BNL). The 
apparatus is described in detail elsewhere [12] and is 
shown schematically in Figure 1. The MPS was upgraded 
for E852 with the inclusion of a cylindrical drift cham- 
ber (TCYL) [13] and a cylindrical cesium iodide (CSI) 
detector [14], both concentric with a LH2 target. Down- 
stream of the target were proportional wire chambers 
(PWCs) for triggering and tracking, a lead/scintillator 
array (DEA) and a 3000-element lead glass electromag- 
netic calorimeter (LCD). The CSI calorimeter consisted 
of 180 segments - 10 longitudinal segments and 18 az- 
imuthal segments. The DEA had the geometry of a pic- 
ture frame with a rectangular central hole. A scintillation 
counter (CPV) was placed just upstream of and matched 
the geometry of the DEA. The lead glass blocks of the 
LGD [15-18] were each 45 cm long with a transverse area 
of 4 X 4 cm^ and were arranged in a 284 x 172 cm^ stack 
with a 2 block x 2 block beam hole to allow for passage 
of non-interacting beam particles. The DEA and CSI de- 
tected low-energy photons that miss the LGD. The LGD 
was used to measure energies and positions of photons 
from tt" and ij mesons. The MPS included an analyzing 
magnet with a 1 T field and drift chambers for tracking. 
The LH2 target was 30 cm long and 6 cm in diameter. 
The beam was an 18.3 GeV/c negatively charged beam 
consisting mainly (95%) of tt" and had a momentum bite 
A.p/p of 3% and a momentum resolution Sp/p oi 1%. The 
target to LGD distance was approximately 5 m. The in- 
teracting beam trigger required no signal from a small 
scintillation counter placed in the path of the beam in 
coincidence with beam defining counters upstream of the 
MPS. Event triggers included all-neutral, all charged and 
a mix of neutral plus charged particle topologies. The 
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TABLE I: Summary of requirements applied to the data set 
leading to tt~p — » 7r~7r°7r°p. 



Cut requirement 


Number of Events 


Satisfy topology trigger" 


123,748,800 


Require four photons" 


13,737,265 


Confidence level > 0.20" 


9,813,524 



Vertex location'' 8,325,895 

Recoil proton angle and momentum'' 6,620,596 

CSI and DEA maximum energies'' 4,939,661 

Photon separation'' 9,579,426 

Final sample (all requirements) 3,025,981 

"Initial cuts, applied in the order shown. 
''These cuts applied separately, after initial cuts. 



TABLE II: Suniniary of rctiuireineiits applied to the data set 
leading to tt ]) — • - tt - ' j>. 



Cut requirement 


Number of Events 


Satisfy topology trigger" 


78,659,511 


No photons" 


16,796,457 


Confidence level > 0.20" 


10,157,455 


Beam hole'' 


8,038,452 


Vertex location'' 


8,038,452 


Recoil proton angle and momentum'' 


7,030,216 


CSI and DEA maximum energies' 


4,822,330 


A"'"^ elimination'' 


9,749,031 


Final sample (all requirements) 


2,585,776 



"Initial cuts, applied in the order shown. 
''These cuts applied separately, after initial cuts. 




FIG. 2: Scatterplot of one di-photon mass against the other 

di-photon mas (three such pairings per event) for events sat- 
isfying the 7r^7r''7r''p final state trigger and with four recon- 
structed photons in the LCD. The inset shows the di-photon 
mass spectrum with a fit to a Gaussian plus a polynomial. 
The n° mass resolution is 10.4 MeV/c^. 



neutral particle requirement depended on energy and/or 
multi-photon effective mass information from the LGD. 

B. Event reconstruction and selection 

We refer to the w~w'^Tr'^p system as the neutral mode 
and the 7r~7r~7r+p system as the charged mode. In Ta- 
bles I and II the selection requirements and resulting 
event sample sizes are summarized for the neutral and 
charged modes respectively. In what follows we discuss 
details of these requirements [19]. 

1. Initial cuts 

The initial requirement was satisfaction of the topol- 
ogy trigger. For both modes a recoil charged track 
was required through the tracking chamber (TCYL) sur- 
rounding the target. The neutral mode required a single 
forward-going charged particle in the PWCs along with 
energy deposition in the LGD corresponding to multi- 
photon effective mass greater than 200 MeV/c^. For 
the charged mode, three forward-going charged particles 
were required in the PWCs. 



Four photons reconstructed in the LGD was the next 

requirement for the neutral mode. The scatterplot of 
one di-photon mass against the other di-photon mass is 
plotted in Figure 2 - there are three such pairings per 
event in this plot. A clear peak corresponding to tt'^tt'^ 
production is observed. The inset of Figure 2 shows a fit 
of the 27 mass spectrum to a Gaussian and a polynomial. 
The 7r° mass resolution is 10.4 MeV/c^. For the charged 
mode further cuts required no photons in the event. 

The distributions in missing mass squared recoiling 
against the 37r system for the neutral and charged modes 
are shown in Figure 3 (unshaded). Although a charged 
track corresponding to the recoil proton is required, its 
momentum is not measured. A kinematic fitting pro- 
gram, SQUAW [20], is used to vary the measured mo- 
menta of the the three pions and incoming beam, within 
errors, to constrain the missing mass to be rrip. For the 
neutral mode, the two di-photon effective masses were 
constrained to the 7r° mass. 

After the kinematic fit was applied the distribution in 
the confidence level [C.L.) of the fit was observed to be 
flat for C.L. > 0.2 (see Figure 4), and this cut was applied 
to both neutral and charged modes. 
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FIG. 3: Distribution in missing mass squared recoiling against 
the Stt system for the (left) 7r~7r'^7r'' mode and the (right) 
7r~7r~7r"'" mode. The unshaded distribution is before any cuts. 
The shaded distribution is for events that survived the con- 
fidence level cuts and all the cuts listed in Tables I and II. 
The missing mass is calculated before the kinematic fitting. 
Results of fits of the shaded distributions to a Gaussian are 
shown. 
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FIG. 5: Distribution in the n'^p effective mass for events in the 
charged mode satisfying the initial cuts. The shaded distri- 
bution shows the events accepted by the m^+p > 1.5 GeV/c^ 
requirement. 
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FIG. 4: Distribution in confidence level after kinematic fitting 
for the (left) ■k~7v'\'^ mode and the (right) 7r~7r~7r^ mode. 
Events were required to have a confidence level greater than 
0.20. 



2. Other requirements 

The PWCs used to define the trigger topology were 
inefficient in the region where the beam intercepted the 
chambers. In order to eliminate uncertainties in the over- 
all acceptance, events were rejected if a charged particle 
was within a 2.5 ct region around the beam centroid at the 
plane of each chamber. The beam centroid and widths 
were determined from the data using measurements from 
beamline elements. Clearly, the fractional loss is greater 
for the charged mode, with three tracks, than the neutral 
mode with a single track. 

The vertex of the reaction was required to be well- 
contained within the target. The z-position (along the 
beam direction) of the vertex was required to be within 
the central 28 cm of the 30-cm long target. In addition, 
the position of the vertex in the transverse plane was 
required to be within 2.5 a of the nominal beam position, 
resulting in an elliptical cut. This vertex cut, along with 
the confidence level cut, eliminated events with Ks — > 



TT+TT or Ks — > tt'^tt'^. 

The momentum and direction of the recoil proton can 
be inferred from the results of the kinematic fit. This 
inferred recoil direction was required to be within 20° of 
the direction measured by the TCYL. Furthermore the 
computed magnitude of momentum of the recoil proton 
track had to be consistent with the recoil proton having 
sufficient momentum to escape the LH2 target. The dis- 
tance the recoil proton travelled through the target was 
estimated based on the vertex information and the angle 
of the proton as measured by TCYL. 

The resolution of the charged particle tracking and the 
LCD were not sufficient to rule out events with low en- 
ergy TT^'s whose decay photons did not enter the LCD. 
Thus the CSI and DEA were used to eliminate such 
events based on energy deposition in these calorimeters. 
Information from the TCYL was used to track charged 
particles into segments of the CSI detector. All segments 
without a charged particle passing through them were 
required to register less than 20 MeV of deposited en- 
ergy. Also, if no charged particle passed through the 
CPV counter in front of the DEA, the DEA was required 
to register no signal. 

The event selection requirements associated with the 
beam hole, vertex location, recoil proton and energy de- 
position in the CSI and DEA were applied identically 
to the neutral and charged modes. The remaining cuts 
applied were particular to the two different modes. For 
the neutral mode it was required that the minimum pho- 
ton cluster separation was greater than 8 cm to insure 
reliable photon reconstruction. An additional cut for 
the charged mode required the Tr+p effective mass be 
greater than 1.5 GeV/c^ to eliminate a small signal due 
to A''"+(1236) TT+p as shown in Figure 5. 

After all cuts are imposed 3.0M events remain in the 
neutral mode and 2.6M events remain in the charged 
mode. The resulting missing mass squared distribution 
(before the kinematic fit) is shown in the shaded distri- 
butions of Figure 3. 
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FIG. 6: (a) Uncorrected distribution in the square of the four- 
momentum transfer, t, from the incoming tt" to the outgoing 
(Stt)^; (b) (37r)^ effective mass distribution; (c) tt^tt" effec- 
tive mass distribution; and (d) 7r°7r° effective mass distribu- 
tion for the neutral mode. 
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FIG. 7: (a) Uncorrected distribution in the square of the four- 
momentum transfer, t, from the incoming tt^ to the outgoing 
(37r)~; (b) (37r)~ effective mass distribution; (c) 7r~7r~ effec- 
tive mass distribution; and (d) 7r~7r"'" effective mass distribu- 
tion for the charged mode. 



C. Effective mass and t distributions 

1. Uncorrected distributions 

Distributions of the square of the four-momentum 
transfer, t, from the incoming tt~ to the outgoing (37r)~ 
system as well as the Stt and 2?: effective mass distribu- 
tions are presented in Figures 6 and 7 for the events pass- 
ing the cuts listed in Tables I and II. These distributions 
have not been corrected for acceptance. The distribution 
in t for the neutral mode is shown in Figure 6(a). The 
partial wave analysis, discussed below, was carried out 
for 12 bins in t defined in Table III. The 7r-7r°7r° effec- 
tive mass distribution is shown in Figure 6(b). A peak 
at the 02(1320) is observed as well as an enhancement in 
the vicinity of the 7r2(1670). The tt^tt" mass distribution 
in Figure 6(c) shows clear indication of the p(770) and 
the 7r"7r" in Figure 6(d) shows the /2(1275). 

The distribution in t for the charged mode is shown 
in Figure 7(a). The 7r~7r~7r+ effective mass distribu- 
tion is shown in Figure 7(b). Peaks at the 02(1320) and 
7r2(1670) are clearly observed. The 7r~7r~ mass distribu- 
tion in Figure 7(c) is featureless and the TT~ir~^ in Fig- 
ure 7(d) shows the p(770) and the /2(1275). 

Figure 8(a) shows the 7r~7r°7r'^ effective mass distribu- 
tion for three regions of t that we label tl, t6 and t8 and 
define in Table III. Figure 8(b) shows the 7r~7r~7r+ ef- 
fective mass distribution for the same three regions of t. 
Clearly the overall shape of the 37r effective mass spec- 
tra strongly depends on t. These mass spectra are not 



TABLE III: Definition of the bins in t - square of the four- 
momentum transferred from the incoming 7r~ to the outgoing 
(37r)~ system - used in this analysis. 



t-Bin 


Range in {GeV/cf 


tl 


0.08 to 0.10 


t2 


0.10 to 0.12 


t3 


0.12 to 0.14 


t4 


0.14 to 0.16 


t5 


0.16 to 0.18 


tG" 


0.18 to 0.23 


tl 


0.23 to 0.28 


t8 


0.28 to 0.33 


t9 


0.33 to 0.38 


tio 


0.38 to 0.43 


tn 


0.43 to 0.48 


tl2 


0.48 to 0.53 


''The results presented 


in this paper are for this t bin unless noted 



otherwise. 



corrected for acceptance. 

The various 27r effective mass spectra also depend on 
37r mass. In Figure 9 we show several 27r spectra with 
the requirement that the 37r mass lie in the a2(1320) mass 
region (1.20 < ms^ < 1.44 GeV/c^) or in the 772(1670) 
region (1.54 < ms^ < 1.80 GeV/c^). The 7r~7r° mass 
spectrum (Figure 9(a)) and tt^tt^ mass spectrum (Fig- 
ure 9(c)) are both dominated by the p{770) when the 37r 
mass is in the 02 region. The tt'^tt'^ mass spectrum (Fig- 
ure 9(b)) and 7r~7r"'" mass spectrum (Figure 9(d)) both 
show strong signals for the /2(1270) when the 37r mass 
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FIG. 8: Uncorrected (Stt)" mass distribution for the (a) neu- 
tral mode and (b) charged mode for three t regions: ti , te and 
ts as defined in Table III. 
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FIG. 9: Uncorrected 2n mass distributions for selected regions 
of the Stt effective mass, (a) n~Tv'^ mass in the a2 region; 
(b) tt^tt" mass in the 712 region; (c) 7r~7r^ mass in the 02 
region; and (d) tt^tt^ mass in the -K2 region. 



is in the 7r2 region with the tt 7r+ mass spectrum also 
showing a p{770) resonance signal. 



2. Mass and t acceptance and resolution 

The dependence of the experimental acceptance on rel- 
evant kinematic variables was estimated by generating 
Monte Carlo events for the two modes. Events were 
generated with an exponentially damped distribution in 
t, based on what is observed in the data. The gener- 
ated distribution in Stt mass was uniform from SttItt to 
2.5 GeV/c^. The 2tt masses were chosen to uniformly 
populate the Stt Dalitz plot at fixed Stt mass and the 
distributions in the relevant decay angles were also gen- 
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FIG. 10: (a) Stt mass acceptance as a function of Stt mass; 
(b) t acceptance as a function of t; (c) 5m - mass resolution as 
a function of 3n mass; and (c) St - t resolution as a function 
of t mass - all for both neutral and charged modes. 



erated uniformly. 

The response of each detector component to a charged 
pion or photon from tt*^ decay was simulated along with 
relevant smearing of momenta or energies. The detector 
responses were used as input to the same software used 
to reconstruct tracks and photons from actual data. The 
output was then passed on to the kinematic fitting soft- 
ware. The cuts summarized in Tables I and II were then 
imposed on the simulated event samples. 

The acceptance as a function of Stt effective mass, 
for both neutral and charged modes, is shown in Fig- 
ure 10(a). The acceptance as a function of t for both 
modes is shown in Figure 10(b). The mass resolution, 
Sm, as a function of Stt effective mass and t resolution, 
dt, as a function of t for both modes is shown in Fig- 
ure 10(c) and 10(d) respectively. 

Figures 11(a) and 11(b) show the comparison between 
the acceptance-corrected tt^tt^tt^ and 7r~7r~7r^ mass dis- 
tributions for two regions in t, tl and t8 (defined in Ta- 
ble III). 



III. PWA METHODOLOGY 

A. Overview 

The production of the (Stt)" system in the reaction 
TT~p — > (37r)~p is described as a coherent and incoherent 




(a) resonance: X decay (b) isobar: R decay 
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FIG. 11: Acceptance corrected -k~-k tv and 7r~7r~7r^ effective 
mass distributions for two regions in t: (a) tl and (b) t8. 
These t regions are defined in Table III. 
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FIG. 12: Partial wave analysis of the 3tt system in the isobar 
model. The state X is characterized by spin (J), parity (P) 
and charge conjugation (C). It decays at point 1 into a di-pion 
resonance i?^^ (the isobar) and a bachelor n. The di-pion has 
spin S. The angular momentum between the isobar and the 
bachelor vr is L. At point 2 the di-pion resonance decays into 
TTa and TTft. More details are given in the text. 



sum of partial wave amplitudes. The production through 
one such partial wave amplitude, assuming the isobar 
model, is shown schematically in Figure 12. The state X 
described by the partial wave decays at point 1 into a di- 
pion resonance Rt^t^ (also referred to as the isobar) and a 
bachelor tt followed by the decay of the di-pion resonance 
at point 2 into ttq and 7rt,. The state X is characterized 
by J^'-'' (SL), where J, P and C are the spin, parity 
and charge conjugation of the Stt system, respectively, M 
is the spin projection along the z axis and e represents 
symmetry (reflectivity) of the Stt system under reflection 
in the production plane. Viewed in terms of some ex- 
change mechanism (as in Figure 12) e = + corresponds 
to natural parity exchange and e — ^ corresponds to 
unnatural parity exchange. 

The decay at point 1 is described in the Gottfried- 
Jackson frame [21] which is the rest frame of the Stt sys- 
tem with the beam direction defining the z-axis and the 
normal to the production plane (specified by the momen- 
tum vectors of the beam and recoil proton) defining the 
y-axis. In this frame the angles of the bachelor pion are 
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FIG. 13: The decay angular structure for a particular partial 



wave where the resonance (X) has J 



PC 



and decays 



into /2(1275)7r where the isobar is the /2(1275) and the an- 
gular momentum between the isobar and bachelor tt is L = 2. 
(a) The correlation between the decay angles {Pgj,4'G.j) of 
the resonance in the Gottfried-Jackson frame and (b) The 
correlation between the decay angles {6H,<t>H) of the isobar 
in the helicity frame. 



denoted by Oqj and (pcj- The angular momentum be- 
tween the bachelor pion and i?7rir is and the spin of 
i?7r7r is "S*. The decay at point 2 {Rtt-k — > T^aT^b) is de- 
scribed in the helicity frame - the rest frame of the TTaTTf, 
system with the boost direction to that frame defining 
the z-axis. In this frame the angles of one of the decay 
pions is denoted by Oh and <j)H- The angular distribu- 
tions for the isobar decay are given by the spin of i?7r7r 
and the line shape as a function of m^^ by a relativistic 
Breit-Wigner function with a Blatt-Weisskopf factor [10] 
with resonance parameters given by the Particle Data 
Group [22]. In this analysis the S'-wave tttt isobars are 
parameterized according to the prescription described in 
reference [10]. 

Figure 13 shows the correlated angular decay struc- 
ture between 0gJt4'gj a-nd OH,(t>H for one particular 
partial wave corresponding to the decay of resonance X 
with J^*-^ = 2 ^ into /2(1275)7r where the isobar is the 
72(1275) and the angular momentum between the isobar 
and bachelor tt is L = 2. A typical PWA involved fit- 
ting coherent or incoherent sums of many such correlated 
distributions to the observed data, while taking into ac- 
count the acceptance of the experiment. The acceptance 
in the relevant decay angles for the Stt effective mass in 
the 7r2(1670) mass region is shown in Figure 14. 

The observed intensity as a function of 7713^ and t is 
written as: 



I{m3^,t,T) = 77(r)^ 



^a|(TO3^,<)A^(T) 



(1) 



where a|(m37r,t) is the production amplitude, AI{t) is 
the decay amplitude and r = {9Gj,<t>GJ,QH,4>H,mTrTr} 
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FIG. 14: The acceptance in decay angles for the 3n effective 
mass in the 7r2(1670) mass region. The angles 6 and (f> are 
measured in the Gottfried- Jackson frame for the resonance 
decay and in the helicity frame for the isobar decay. Accep- 
tances are shown for (a) cos 6 and (b) <f) for both the charged 
and neutral modes of the Stt system. 



is the set of kinematic variables describing the decay of 
the resonance and isobar. The index b enumerates the 
spin variables of the partial wave and e the reflectivity, 
as defined below. The AI{t) involve products of I^ma 
functions for the two decays and a function describing 
the line shape of the isobar. In detail: 



i-JML 



ir) 



V(2i + 1)(25+1) 

X 

x{LOSX\J\) 

xFL(px)Fs(pR)BWim^^) + (1^2) 

(2) 

where 1,2 refer to the two identical pions (tt^'s in the 
charged mode and 7r°'s in the neutral), J,M are the 
spin and projection of the resonance, S is the spin of the 
isobar and L is the orbital angular momentum between 
the isobar and the bachelor pion. The functions Fl{px) 
and Fs{pr) are the barrier factors as a function of the 
breakup momentum, px and p^, for X and the isobar 
respectively. BW{mTr,r) is the Breit Wigner function for 
the isobar. In the reflectivity basis: 



A'jML = e(M) [AjML - eP{-iy-^Aj_ML] 



(3) 
1/2, 



where P is the parity of the resonance, 0(0) 
e(M > 0) = ^yTj2 and e(M < 0) = 0. 

The goal of the PWA is to extract the production am- 
plitudes by fitting the correlated decay angular distribu- 
tions to the intensity described by equation 1 in bins of 
tostt and t. In order to account for detector acceptance, 
Monte Carlo simulations were used to find the acceptance 
function r]{T) in terms of the set of kinematic variables 
T. Thus the data were fitted to the decay amplitudes 
modified by the acceptance function. 

The PWA formalism used here is identical to that used 
in [10] and [23] and is also briefly described in inter- 
nal notes [19]. The PWA software was developed at 



Brookhaven Lab and Indiana University (lU) [19]. The 
lU software was optimized for running on a 200-processor 
computer cluster (AVIDD) [24] allowing systematic stud- 
ies involving many PWA fits, varying isobar parameters 
and using different wave sets. The lU programs were also 
used to analyze the data presented in [10] and produced 
consistent results (see Appendix A for a list of technical 
notes available online). 



B. Selecting a partial wave set 

1. General criteria 

As noted earlier, it is important to include in the analy- 
sis a set of partial waves sufficient to describe the physics 
and to have well understood criteria for selecting this set. 
One technique, following [8], is to start with some par- 
ent set of waves and examining the change in likelihood 
resulting from sequentially removing waves. 

Another arbiter of wave set sufficiency is the com- 
parison of moments, H{imm'), of the Df^^,{Q) func- 
tions as calctilated directly from the data and as calcu- 
lated using the results of PWA fits. A PWA using a 
sufficient set of partial waves will reproduce all angu- 
lar distributions of the data as evidenced by the agree- 
ment of the observed and calculated moments. We define 
H{£mm') = J I{fl)Df^^,{fl)dfl where fl represents the 
Euler angles of the Stt system and the intensity /(fi) is 
determined directly from experiment or computed using 
the results of the PWA fits. Both of these techniques 
were used in this analysis. 



2. Defining a parent set of waves 

The TTTT isobars included in this analysis are the a, 
/o(980), p(770), /2(1275) and P3(1690) where a is meant 
to indicate a 5- wave tttt system as described in [10]. 
In addition, a backgrotmd wave is inchided. The back- 
ground wave is characterized by a uniform distribution in 
the relevant decay angles and is added incoherently with 
the other waves. We restrict the parent wave set to those 
with J < 4, M < 1 and 5* < 3. 



3. Removing waves from the parent set 

Waves were sequentially removed from the parent set 
and the resulting change in likelihood (jC) examined [19]. 
Three passes were used. First, waves with an L=0 iso- 
bar (tttt 5- wave and /o) were examined. Subsequent 
passes examined L = 2 (/2) and L = 1 (p) isobar par- 
tial waves. A PWA with the wave in qticstion removed 
was performed and the resulting likelihood compared to 
PWA including the wave. If the change in ln£ was less 
than 40 the wave was removed. Figure 15 shows the 
change in InC when a significant, but not a dominant. 
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TABLE IV: Partial waves under consideration in this analysis. 
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FIG. 15: Partial wave intensity (filled circles) as a function 
of Stt mass (lower plot) shown along with the change in In jC 
(shown above the intensity plot - read the right axis scale) 
that results when this partial wave is removed from the parent 
set of partial waves, (a) 2 ^1^ f2nD wave - a significant wave 
and (b) 2 f2nD wave - an insignificant wave. 



the low wave set alone [L] or in both sets [B]. 



wave (2 +l+/27rZ? wave) is removed compared to when 
an insignificant wave (2 ^l~/27r D wave) is removed. 

Certain partial waves that would have been removed 
by this selection criterion were kept because the exis- 
tence of signals in these waves is under consideration. 
The negative reflectivity exotic 1 partial waves were 
kept even though they would have been chosen for re- 
moval by the above criterion. The positive reflectivity 
exotic partial wave failed the selection criterion every- 
where above 1.4 GeV/c^ and exhibited large fluctuations 
in intensity below this mass. This partial wave was also 
kept and its contribution to the likelihood is further ex- 
amined below using the final set of partial waves. Finally, 
the 2++0^/97r i'-wave (which failed the selection test) was 



kept to serve as a possible reference wave for evaluation 

of phase differences of negative reflectivity partial waves. 
Based on these criteria the compiled set includes 35 waves 
and a background wave. We refer to this as the high-wave 
set. 

The analysis reported in reference [10] used a wave set 
consisting of 20 waves including the background wave - 
we refer to this as the low-wave set. Table IV lists the 
waves used in low-wave and high-wave sets. The waves 
included only in the low- wave set are labeled with [L] in 
Table IV while the waves included only in the high-wave 
set are labeled with [H] and those included in both sets 
are labeled with [B]. 
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FIG. 16: Comparison of the (a) H(201) and (b) H(420) mo- 
ments as computed directly from data and from PWA fits for 
the low and high wave sets for the 7r~7r''7r'' channel. In (c) 
the A'^ (differences squared divided by errors squared summed 
over all mass bins and divided by the number of mass bins) 
is shown for various moments. 



FIG. 17: Comparison of the (a) H(201) and (b) H(420) mo- 
ments as computed directly from data and from PWA fits for 
the low and high wave sets for the tt^tt^tt^ channel. In (c) 
the (differences squared divided by errors squared summed 
over all mass bins and divided by the number of mass bins) 
is shown for various moments. 



4- Using moments 



The PWA solutions obtained with the two wave sets 
were used to compute moments that were then compared 
directly with data. In Figure 16 we show the iJ(201) and 
7f (420) moment comparisons as a function of Stt mass for 
the 7r~7r'^7r° mode for fow-wave and high- wave set PWA 
fits. We also show the difference between data and PWA 
calculations of the moments as (summing differences 
squared divided by errors squared summed over all mass 
bins and divided by the number of mass bins) for various 
moments. Similar plots are shown in Figure 17 for the 
7r~7r~7r+ mode. The moments calculated using the PWA 
results from the high wave set have better agreement with 
experimental moments, but there are moments (such as 
H{201) - shown in part (a) of Figure 17) for which agree- 
ment is not achieved. This may be due to the inherent 
inadequacy of the isobar model in describing the under- 
lying production mechanism. 

An examination of Figures 16(c) and 17(c) show that 
for all the moments, the moments calculated using the 
high-wave set have better agreement with data than do 
moments calculated using the low-wave set. 



IV. PWA RESULTS 

The PWA results presented in this section were for 
data with 0.18 < |t| < 0.23 (GeV/c)^ referred to as t6 in 
Table III. The analysis was also carried out in 11 other 
bins in \t\ - some results of which will also be summarized 
below. 

Figure 18 shows the sum of all waves in the high-wave 
set along with the background wave for both the charged 
and neutral modes. The fitting procedure constrains 
the sum of all the waves to agree with the acceptance- 
corrected Stt mass distribution. The background wave 
(uniform in all angles) is added incoherently with the 
other waves. As can be seen, the contribution from the 
background wave is small (of order 1%). 

There is a systematic disagreement in the relative nor- 
malization between the two modes resulting from an over- 
all per track inefficiency of 5% in the trigger PWC's. 
The overall inefficiency has no bearing on the features of 
the acceptances and was thus not incorporated into the 
Monte Carlo corrections, but results in the neutral mode 
having a yield consistently 25% higher than the charged 
mode. The neutral mode required only one forward-going 
charged particle in the trigger while the charged mode re- 
quired three. 

In the following we discuss intensities and phases of 
individual partial waves. A partial wave intensity corre- 
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FIG. 18: The sum of a partial waves and the background wave 
for both the charged and neutral modes. 



spends to a single term of equation (1) which is a product 
of the fitted partial wave amplitude, a^(m37r,t), and the 
isobar amplitude Ajml{t) ■ Similar to the procedure fol- 
lowed in references [9, 10], the intensities are integrated 
over acceptance corrected data, and are normalized to 
the total number of events. 

The coherent sum of the 1^"*" /jtt- waves is shown in 
Figure 19 for the 7r^7r^7r+ and tt^vt^tt^ modes. For each 
the results from the high wave set and low wave set are 
shown. The 1+"*" wave is a dominant wave and the results 
of the PWA for the two wave sets (high and low) are in 
good agreement. The yield for the two charged modes is 
also consistent with expectations from isospin for the p7r 
system. 

The 1"'"'" pTT D-wave is not included in the high- wave 
set, since it did not meet criteria for inclusion, but it is 
included in the low-wave set. The observed ratio of 1++ 
/97r S-wave to 25-wave in the low-wave set is consistent 
with that of reference [10]. 

Figure 20 shows the 2++1+ D-wa.ve pn partial wave 
as a function of ms^r for both the neutral and charged 
modes. A clear 02 (1320) signal is observed. For a pure 
pn isovector resonance, we expect an equal number of 
charged mode and neutral mode events. This wave is 
clean and robust and will be used for interferometry to 
study the phase motion of other waves relative to this 
wave. 

Figures 21 and 22 show the intensities for the 2 '"0+ S- 
wave /27r and 2 ^0^ D-wave /27r waves respectively for 
both the charged and neutral modes. Note that the rela- 
tive yields between charged and neutral modes, for these 
/27r decays are consistent with expectations from isospin, 
after taking into account the overall 25% reduction in 
yield for the charged mode. These intensities are also 
compared with those measured by Daum et. al. [8] in 
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FIG. 19: Coherent sum of the 1"*"^ p7r-waves for the (upper) 
7r~7r~7r+ and (lower) ■k'tt'^-k" modes. For each the results 
from the high wave set and low wave set are shown. 



a similar amplitude analysis of diffractive Stt production 
in Tr~p interactions at 63 and 94 GeV/c^. The Daum 
results have been arbitrarily scaled by a factor of 5. 

The S and D fi-n modes both show a resonant-like 
structure but not at the same peak position. The S wave 
peaks at approximately 1.6 GeV/c^ while the D wave 
peaks at approximately 1.8 GeV/c^. This shift in the 
peak position has been observed in other analyses [8, 10]. 
This has been interpreted as the /27r decay mode of two 
different resonances. The PDG places the 7r2(1670) 2 ^ 
resonance at a mass of 1.67 GeV/c^ based on fits to the 
S wave mode alone. 

Figure 23 shows the 4++0+ G-wave p7r intensity as a 
function of ma^n- for both the neutral and charged modes. 
Comparing Figures 23 and 20 the yield of 04(2040) is 
about 3% that of the 02(1320) indicating the sensitivity 
of the PWA within the isobar model for finding states 
with relatively low production cross-sections. 
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THE EXOTIC 1 + WAVE 



B. Leakage of the 7r2(1670) and the instability of the 
exotic wave 



A. Comparing the high and low wave sets 

The 1 ''1+ P-wave pn partial wave and its interfer- 
ence with the 2++ wave are shown in parts (a) and (b) 
of Figures 24 and 25 for the neutral mode (7r~7r°7r°) and 
charged mode (7r~7r~7r+). When the low- wave set is used 
in the fit, an enhancement is observed for the exotic 1 ^ 
wave in the Sir mass region around 1.6 GeV/c^, consis- 
tent with the observation of [9, 10]. The enhancement 
disappears when the high wave set is used. This effect 
is observed for both Stt channels. The phase of the ex- 
otic wave relative to the dominant 2++ wave (see Fig- 
ures 24(b) and 25(b)) is similar for both wave sets and 
dominated by the phase of the 02 (1320). The apparent 
phase motion beginning around 1700 MeV/c^ is likely 
due to the 02(1700) [22], but the unreliability of the ex- 
otic intensity makes this an inappropriate place for its 
study. 



One source of leakage into the exotic wave for the 
PWA carried out with the low- wave set is the 7r2(1670). 
In Figure 26 the coherent sum of the intensities for the 
2 '"0"'" F-wave pn, 2 '"1"'" F-wave pn and 2 '"1+ P-wave 
yOTT waves for the charged and neutral modes are shown. 
These waves are allowed decay modes of the 7r2(1670). A 
clear peak is seen at the 7r2(1670) mass, the same mass 
region where the purported exotic meson was seen using 
the low- wave set. The three waves shown in Figure 26 
were not included in the low- wave set. 

The PWA fit was repeated for the high wave set but 
with the 2 ^0+ P-wave pn, 2 ""1+ P-wave pn and 
2 '"1+ P-wave pir waves removed. The resulting 1 ^ P- 
wave pTT intensities for the neutral mode for this modified 
high wave set fit, along with the unmodified high wave 
set and low wave set, are shown in Figures 27. The corre- 
sponding results for the charged mode are shown in Fig- 
ure 28. The agreement between the modified high wave 
set and the low wave set shows that the 2 P-wave pir 



13 



> 
(J 

in 
rN 
o 
d 



^ 1000- 



2 +0+ D-wave i^K 

I I I I L_ 



• high wave 
low wave 
D Daum 




' — 1 1 r 

1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 
J I I I I L 



> 

in 

o 
p 



2 1000- 




M[37i] GeV/c^ 

FIG. 22: 2 ^0"*" D-wave /27r partial wave intensity for the 
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FIG. 24: (a) The l""*"!^ P-wave partial wave in the neu- 
tral mode (7r~7r''7r'') for the high-wave set PWA and the low- 
wave set PWA and (b) the phase difference A$ between the 
2^"*" and 1 ^ for the two wave sets. 
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FIG. 25: (a) The 1 ^1"*" P-wave pTv partial wave in the 
charged mode (7r~7r~7r"'') for the high-wave set PWA and the 
low-wave set PWA and (b) the phase difference A$ between 
the 2+"*" and 1""'' for the two wave sets. 
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FIG. 26: Sum of the intensities for the 2 ^0 P-wave pn, 
2~ + i+ P-wave pvr and 2 '"1^ P-wave pn waves for the 
charged and neutral modes. These waves were included in 
the high- wave set but not in the low- wave set. 
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FIG. 27: The positive refieetivity 1^^ P-wave pTv intensity 
for the neutral mode for the high wave set (filled circles) , the 
modified high wave set (filled squares) and the low wave set 
(open circles). In the modified high wave set the two 2 F- 
wave pTT waves and the 2 P-wave piT wave were removed. 
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FIG. 29: The 2"'""'" intensity for charged Stt mode for the high 
wave set (filled circles) and low wave set (open circles). The 
4+"'" intensity (filled triangles) for the high wave set is also 
shown. 
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the high wave set clearly shows the 04(2040). Figure 29 
shows a semi-log plot of the the 2++ intensity for charged 
Stt mode for the high wave set (filled circles) and low wave 
set (open circles). The 4++ intensity for the high wave set 
is also shown (filled triangles). The 02(1320) is dominant 
in the 2+"'" intensity in both high and low wave sets. The 
2++ intensity for the low wave set shows a resonance- 
like structure near 2 GcV/c^ in low wave set but not the 
high wave set. However the 2+"'" low wave set intensity 
agrees well with the 4++ high wave set intensity. In this 
case, leaving out the 4++ in the low wave set forces a 
enhancement at 2 GeV/c^ in the 2++ wave. 



FIG. 28: The positive reflectivity 1 P-wave pn intensity 
for the charged mode for the high wave set (filled circles) , the 
modified high wave set (filled squares) and the low wave set 
(open circles) . In the modified high wave set the two 2 F- 
wave piT waves and the 2 l""" P-wave pir wave were removed. 



D. Further discussion of the significance of the 
exotic wave 



decays and the 2~+l+ P-wave pn decay of the 772 are the 
source of leakage leading to the observed exotic peak in 
the low wave set. 



C. Evidence for leakage of the O4(2040) 

A similar phenomenon is observed in the 2++ wave 

in the low wave set where the 4"'"'"0"'" G-wavc pn is not 
included. As discussed above, the 4++0+ G-wave p7r in 



The contribution to the In C of the positive reflectivity 
exotic partial wave was examined by removing this wave 
from the high wave set and comparing the resulting likeli- 
hood with that of the original high wave PWA. Figure 30 
shows the resulting change in likelihood. The sign con- 
vention of this difference is such that positive (negative) 
values indicate inclusion of this partial wave improved 
(degraded) the fit quality. As can be seen from figure 30 
there is no region of 37r effective mass where this partial 
wave is consistently required by the fit. This behavior is 
seen in all t-bins presented here. 
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FIG. 30: The change in InL resuhing from the removal of 
the 1 '"l"'"(p7r)P-wave partial wave. See the text for the sign 
convention and discussion. 
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FIG. 32: Variation of the 2 '"0"'" 5-wave /27r intensity for the 
charged mode for the 12 t-bins of Table III, increasing in \t\ 
from left to right and top to bottom. 



shown in Tabic III for the high wave fit in the charged 
Stt. As noted in Table III, the width of the |t| bin for 
the first five bins is 0.02 (GcV/c)^ and 0.05 (GeV/c)^ for 
the other bins. Figures 31, 32 and 33 show the variation 
of the 1 '"1+ P-wave pw intensity, the 2 '"0+ S-wave 
/27T intensity and the phase difference between the 1 ^ 
and 2 amplitudes for the charged 37r mode for the 12 
f-bins. The 2 '"0+ 5-wave /27r intensity is stable across 
the 12 t bins. We note that although the 1~+1+ P- 
wave pn intensity may indicate an enhancement in the 
vicinity of ~ 1.6 GeV/c^ for some of the higher t bins 
in Figure 31, the corresponding phase difference shown 
in Figure 33 is not consistent with resonant behavior. 
Indeed the 2~+0+ S'-wave /27r and 1 '"1+ P-wave pw 
are phase-locked in the 1.6 GeV/c^ region. One possible 
explanation is production of an exotic resonance with 
the same mass, width and phase motion as the 7r2(1670). 
Another is some residual leakage of the 712 (1670) in the 
exotic wave in the high- wave set. 



E. t dependence studies 

1. The exotic wave 

We examined the stability of the exotic wave inten- 
sity and the intensity of other waves across the 12 t-bins 



2. The 2++ and 4++ waves 

The 2++1"'" D-wave pn intensity (shown in Figure 34) 
and the 4++0+ G-wave pw intensity (shown in Figure 35) 

are shown for the 12 t bins. For completeness, in Fig- 
ure 36, we compare the t dependence of the 2++1+ D- 
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FIG. 33: Variation of the phase difference between the 1 
and 2 amplitudes, whose intensities are shown in Figures 31 
and 32 respectively, for the charged mode for the 12 t-bins of 
Table III, increasing in \t\ from left to right and top to bottom. 




M[37C] GeV/c" 



M[37t] GeV/c^ 



M[3jt] GeV/c^ 



FIG. 34: Variation of the 2++1"'" D-waye pn intensity for the 
charged mode for the 12 t-bins of Table III, increasing in |t| 
from left to right and top to bottom. 



wave pir intensity, dominated by the 02 (1320), in both 
Stt modes in this experiment, with other measurements 
of TT-p a2 (1320)p at 18.8 GeV/c[25]. 



VI. STUDIES OF SYSTEMATICS 

A. Effects of varying data selection cuts 

Studies were carried out to determine the effect on the 
PWA of varying the parameters that define the selec- 
tion criteria as described in section II B. Different crite- 
ria were found to effect the results in different ways but 
no significant change in our conclusions was observed. 
The effects of these perturbations on the analysis are de- 
scribed individually below. 

The beam hole cut is nominally 2.5a. Values of 1.5a 
and 3.5a systematically increased and decreased the in- 
tensity of the l+^pTT S-wa.ve below 1.2 GeV/c^ but left 
the mass region above 1.2 GeV/c^ indistinguishable from 
the unperturbed analysis. The changes observed were ap- 
proximately 5% of the intensity. Conclusions regarding 
the low mass behavior of this partial wave should take 
into account this systematic effect. Other partial waves 
were unaffected. 



The cut to remove the A++ in the charged sample 
affected the data only above 2.0 GeV/c^. Removing 
this cut from the analysis increased the intensity of the 
2 ^0+(/27r) 5*- wave by approximately twice its error be- 
tween 2.0 and 2.4 GeV/c^. The systematic effect on other 
partial waves is insignificant when compared with their 
statistical errors. 

The effect of removing the DEA cut was similar to 
using a less restrictive value for the beam hole cut. Small, 
but statistically significant, increases in the intensity of 
the 2 ^0+(/27r) S'-wave were observed between 1.2 and 
1.6 GcV/c2.' 

Increasing (decreasing) the confidence level cut to 0.3 
(0.1) changed the observed intensities everywhere by a 
mass independent scale factor. The same effect was ob- 
served for the CSI cut. Removing this cut increased all 
intensities by a factor of 1.66 but left their shapes un- 
changed. 

In summary, we have investigated the effect of perturb- 
ing the selection criteria summarized in Tables 1 and II 
on the results of the partial wave fit and have shown that 
our results do not depend on the specific values used in 
the criteria. 
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FIG. 35: Variation of the 4"''^0^ G-wave pvr intensity for the 
charged mode for the 12 t-bins of Table III, increasing in \t\ 
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B. 5- wave parametrization 

We used the isobar model, following the same method- 
ology employed in references [9, 10] . This model uses two- 
body amplitudes for the tttt system. The 5 = 1 isovector 
and S — 2 isoscalar amplitudes were parametrized us- 
ing Breit-Wigner amplitudes. There are various = 0, 
isoscalar parametrizations available in the literature. In 
general the tttt amplitude is dominated by a broad a me- 
son, the /o(980) and several resonances above 1 GeV/c^, 
depending on the specific parametrization. The results 
presented in this paper were all based on the parametriza- 
tion used in references [9, 10] . This is based on the 
model of [26] with the possibility of direct /o(980) pro- 
duction; the /o(980) was subtracted from the elastic tttt 
amplitude and added as a Breit-Wigner resonance. We 
have also studied other parametrizations, e.g. based on 
(1) N/D dispersion relations, (2) explicit inclusions of 
the KK channel and/or (3) including direct production 
of scalar resonances above 1 GeV/c^ [27, 28]. All give 
qualitatively identical results. 



VII. CONCLUSIONS 

A partial wave analysis of the tt^tt^tt^ and tt^tt'^tt'^ 
systems has been presented using a data set whose statis- 
tics exceeds published analyses of the Stt systems. The 
ai(1260), 02(1320), 7r2(1670) and a4(2040) mesons were 
observed. The data provide no evidence for an exotic 
'+ meson as had been reported earlier using a 
smaller sample of the Tr^Tr^Tr"*" mode alone [9, 10]. This 
analysis included the study of criteria for defining a set 
of partial waves sufficient to describe the data includ- 
ing a study of the effect on the fit quality by removing 
individual waves and also a comparison of observed mo- 
ments of the data with moments computed from PWA 
solutions. The wave set thus established contained more 
waves than were used in the analysis reporting evidence 
for the exotic meson. Our studies indicate that leaving 
out partial waves corresponding to decay modes of the 
772(1670) lead to a false enhancement in the exotic 1 
wave. This analysis and the analysis of references [9, 10] 
were based on the isobar model. The isobar model is 
known to suffer some limitations, as has been pointed 
by others [30? ]. The so-called Deck mechanism [30? ] 
is particularly relevant to the production of the Stt sys- 
tem and is being applied to these data. Results will be 
presented in a forthcoming paper. 
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FIG. 36: t dependence of the 02(1320) as measured in this 
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APPENDIX A: TECHNICAL NOTES AND 
ONLINE PWA RESULTS 

The authors of this paper have made available techni- 
cal notes and PWA results online [19]. The notes include 
(1) details of the data selection; (2) a comparison of the 



software used in this analysis with that used in the anal- 
ysis of reference [10]; (3) the Monte Carlo details for this 
analysis; (4) a full description of the PWA formalism; 
(5) the moments method technique and definition of an- 
gles; (6) a study of background contributions to partial 
wave intensities; (7) stability of PWA results against 
variations in data selection criteria; (8) acceptance and 
resolution; (9) the determination of the minimal partial 
wave set and (10) S'-wavc paramctrization. In addition, 
PWA results (amplitudes and phases) are available for 
12 bins in t and for various wave sets. 
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